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ABSTRACT 

We present new Spitzer, UKIRT and MMT observations of the blue compact 
dwarf galaxy (BCD) Mrk996, with an oxygen abundance of 12+log(0/H) = 8.0. 
This galaxy has the peculiarity of possessing an extraordinarily dense nuclear 
star-forming region, with a central density of ~ 10 6 cm -3 . The nuclear region of 
Mrk 996 is characterized by several unusual properties: a very red color J — K = 
1.8, broad and narrow emission-line components, and ionizing radiation as hard 
as 54.9 eV, as implied by the presence of the [Oiv] 25.89 /xm line. The nucleus 
is located within an exponential disk with colors consistent with a single stellar 
population of age > 1 Gyr. The infrared morphology of Mrk 996 changes with 
wavelength. IRAC 4.5 /xm images show extended stellar photospheric emission 
from the body of the galaxy, and an extremely red nuclear point source, indicative 
of hot dust; IRAC 8 /xm images show extended PAH emission from the surround- 
ing ISM and a bright nucleus; MIPS 24 and 70 images consist of bright point 
sources associated with the warm nuclear dust; and 160 /xm images map the 
cooler extended dust associated with older stellar populations. The IRS spec- 
trum shows strong narrow Polycyclic Aromatic Hydrocarbon (PAH) emission, 
with narrow line widths and equivalent widths that are high for the metallic- 
ity of Mrk 996. Gaseous nebular fine-structure lines are also seen. A CLOUDY 
model which accounts for both the optical and mid-infrared (MIR) lines requires 
that they originate in two distinct H II regions: a very dense H II region of radius 
~ 580 pc with densities declining from ~ 10 6 at the center to a few hundreds 
cm -3 at the outer radius, where most of the optical lines arise; and a H II region 
with a density of ~ 300 cm -3 that is hidden in the optical but seen in the MIR. 
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We suggest that the infrared lines arise mainly in the optically obscured H II 
region while they are strongly suppressed by collisional deexcitation in the op- 
tically visible one. The hard ionizing radiation needed to account for the [O iv] 
25.89 /iin line is most likely due to fast radiative shocks propagating in an inter- 
stellar medium. A hidden population of Wolf-Rayet stars of type WNE-w or a 
hidden AGN as sources of hard ionizing radiation are less likely possibilities. 

Subject headings: Galaxies: individual: Mrk996; Galaxies: compact; Galaxies: 
dwarf Galaxies: starburst; (ISM:) dust, extinction; ISM: lines and bands; In- 
frared: galaxies 

1. Introduction 

Among blue compact dwarf galaxies (BCDs), the dwarf emission-line galaxy Mrk 996 
(M B = —16.9) occupies a place apart because of the extreme electron density at the center of 
its star-forming region. HST V and / images show that the bulk of the star formation occurs 
in a compact roughly circular high surface brightness nuclear region of radius ~ 340 pc, with 
evident dust patches to the north of it (Thuan et al. 1996). The nucleus (n) is located within 
an elliptically (E) shaped low surface brightness (LSB) component, so that Mrk 996 belongs 
to the relatively rare class of nE BCDs (Loose & Thuan 1985). The extended envelope shows 
a distinct asymmetry, being more extended to the northeast side than to the southwest side, 
perhaps the sign of a past merger. This asymmetry is also seen in the spatial distribution 
of the globular clusters around Mrk 996, these being seen mainly to the south of the galaxy. 

Mrk 996 has a heliocentric radial velocity of 1622 km s -1 which gives it a distance of 21.7 
Mpc, adopting a Hubble constant of 75 km s _1 Mpc -2 and including a very small correction 
for the Virgocentric flow (Thuan et al. 1999). At the adopted distance, 1 arcsec corresponds 
to a linear size of 105 pc. Thuan et al. (1996) found the extended LSB component to possess 
an exponential disk structure with a small scale length of 0.42 kpc. While Mrk 996 does 
not show an obvious spiral structure in the disk, there is a spiral-like pattern in the nuclear 
star-forming region, which is no larger than 160 pc in radius. 

The UV and optical spectra of the nuclear star-forming region of Mrk 996 (Izotov et 
al. 1992; Thuan et al. 1996) show remarkable features, suggesting very unusual physical 
conditions. The He I line intensities are some 2-4 times larger than those in normal BCDs. 
In the UV range, the N ill] A1750 and C ill] A1909 are particularly intense. Moreover, the 
line width depends on the degree of ionization of the ion. Thus low-ionization emission lines 
such as + , S + and N + have narrow widths, similar to those in other H II regions, while 
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high-ionization emission lines such as the helium lines, the ++ and Ne ++ nebular lines and 
all auroral lines, show very broad line widths, > 500 km s -1 . Such correlations of line widths 
with the degree of excitation suggest different ionization zones with very distinct kinematical 
properties. Thuan et al. (1996) found that the usual one-zone low-density ionization-bounded 
H II region model cannot be applied to the nuclear star-forming region of Mrk 996 without 
leading to unrealistic helium and heavy-element abundances. Instead, they showed that a 
two- zone density-bounded H II region model, including an inner compact region with a central 
density of ~ 10 6 cm -3 , some 4 orders of magnitude greater than the densities of normal H 
II regions, together with an outer region with a lower density of ~ 450 cm~ 3 , comparable 
to those of other H II regions, is needed to account for the observed line intensities. The 
large density gradient is probably caused by a mass outflow driven by the large population 
of Wolf-Rayet stars present in the galaxy. The gas outflow motions may account for the 
much broader line widths of the high-ionization lines originating in the dense inner region 
as compared to the low-ionization lines which originate in the less dense outer region. As 
for the high N ill] A 1750, C ill] A 1909 and He I line intensities, they can be understood by 
collisional excitation of these lines in the high density region. In the context of this model, 
the oxygen abundance of Mrk 996 is 12+log(0/H) = 8.0. Adopting 12+log(0/H) = 8.65 for 
the Sun (Asplund et al. 2005), then Mrk 996 has a metallicity of 0.22 solar. Mrk 996 shows 
enhanced helium and nitrogen abundances, which can be accounted for by local pollution 
from Wolf-Rayet stars. 

We present in this paper Spitzer (Werner et al. 2004) mid-infrared (MIR) observations of 
Mrk 996. The extraordinary UV and optical properties of this galaxy, along with the evident 
presence of dust patches in the star-forming region of Mrk 996 on HST optical images, make 
it a prime candidate for our Cycle 1 (PID 3139: P.I. Thuan) Spitzer observations. Our 
entire program consists of spectroscopic, photometric and imaging observations of 23 BCDs 
with metallicities ranging from 1/20 to 1/2 that of the Sun, and its main aim is to study 
star formation in metal-poor environments and to understand how star formation and dust 
properties change as a function of metallicity and other physical parameters. This paper is 
the second of our Spitzer series, the first paper being on Haro 3, the most metal-rich BCD 
in our sample (Hunt et al. 2006). 

We present in §2 our Spitzer IRAC, MIPS and IRS observations of Mrk 996 and their 
reduction. We discuss in §3 new complementary UKIRT near-infrared (NIR) imaging obser- 
vations and optical MMT spectroscopic observations. In §4, we discuss our imaging results: 
the IR morphology of the disk of old stars, and the nature of the very red, bright and dense 
nuclear IR source in Mrk 996. We also discuss the extended PAH emission. In §5, we present 
our spectroscopic results: the PAH features and the IR fine-structure lines. In §6, we use the 
CLOUDY photoionization code (Ferland 1996; Ferland et al. 1998) to model the observed 
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optical and IR emission-line intensities. We show that it is necessary to postulate two H II 
regions: one which is optically visible where the optical lines arise, and one which is opti- 
cally hidden where the main part of the IR line emission arises. We then discuss possible 
sources of hard ionizing radiation - fast shocks, WNE-w stars or an AGN - to account for 
the presence of the [O iv] A25.9yum, line. We summarize our conclusions in §7. 

2. Spitzer Observations 

In the context of our guest-observer program, we acquired Spitzer observations of 
Mrk996: IRAC images (Fazio et al. 2004) at 4.5 and 8 /im, MIPS photometry (Rieke et 
al. 2004) at 24, 70 and 160 /im and low- and high-resolution IRS spectra (Houck et al. 2004). 
The data were obtained during the period from 16 December, 2004 to 3 January, 2005. As 
the data reduction procedures have been described in detail in Paper I (Hunt et al. 2006), 
we give here only a summary. 

2.1. IRAC Imaging 

We acquired four IRAC frames, with small-scale dithering in a cycling pattern, giving 
a total of 120s spent in each of the 4.5 and 8.0 /zm channels. The individual bed frames 
were processed with the Sll.0.1 version of the SSC pipeline which removes the effects of 
dark current, detector nonlinearity, flat field, multiplexer bleeding, and cosmic rays, and 
performs flux calibration. The bed frames were then coadded using MOPEX, the image 
mosaicing and source-extraction package provided by the SSC (Makovoz & Marleau 2005). 
The frames were corrected for geometrical distortion and projected onto a fiducial coordinate 
system with pixel sizes of 1'.'20, roughly equivalent to the original pixels. Standard linear 
interpolation was used for the mosaics. The final coadded 4.5 and 8.0 /xm images are shown 
respectively in the left and right panels of Fig. 1. Fig. 2 shows 4.5 /zm contours superimposed 
on the H5T/WFPC2 F791W optical image of Thuan et al. (1996) (left panel) and on the 
HST 569W - F791W color image from the same authors (right panel). In Fig. 3, we show 
the 4.5/8.0/zm flux ratio image on which have been overlayed the surface brightness level 
contours of the A.bfxm image. To construct the flux ratio image, the two IRAC images have 
been aligned by cross-correlation after sky subtraction. 

We have performed aperture photometry on the IRAC images with the IRAF 1 pho- 
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tometry package APPHOT, taking care to convert the MJy/sr flux units of the images to 
integrated flux. The background level was determined by averaging several adjacent empty 
sky regions. We give the growth curve for both IRAC images (filled circles for the 4.5/im 
image and filled squares for the 8.0/im image) in the left panel of Fig. 4. We have also made 
photometric measurements of the instrument point response functions (PRFs), provided in 
the MOPEX package; the dashed lines show the growth curve expected from a point source, 
normalized to the total flux indicated by horizontal dotted lines in Fig. 4. The 4.5 /im and 
8.0 /im total fluxes are listed in Table 1. 

We have derived radial surface brightness profiles by fitting ellipses to both IRAC images. 
The ellipse centers were held fixed, but the ellipse shape and orientation was allowed to vary 
with radius. The surface brightnesses were put on the Vega magnitude scale by using the 
photometric calibration of Reach et al. (2005). Figure 5 (left panel) shows the radial profiles 
in units of [4.5] and [8.0] magnitudes. The average disk color of 1.5 corresponds to a flux 
ratio 4.5/8.0 of 0.7. The red circumnuclear maximum in the lower panel at a radius of ~ 
4" is an artifact due to the different IRAC point-response functions at the two wavelengths. 
This effect can also be seen in Fig. 3 where the crowding of the 4.5 /im contours from the 
point source coincides with a relatively low 4.5/8 flux ratio. 



2.2. MIPS Imaging 

We have obtained for Mrk 996 a total of 24 frames at 24/im, 56 frames at 70/xm, and 104 
frames at 160/iin, using ramp times of 3s, 10s, and 10s respectively, with 1, 2, and 6 cycles 
in the three channels. The individual bed frames were processed by the Sll.4.0 version of 
the SSC pipeline, which includes dark current subtraction, flat-fielding, and flux calibration. 
As for the IRAC images, we processed the dithered bed frames in the spatial domain with 
MOPEX. Geometrical distortion was corrected before projecting the frames onto a fiducial 
coordinate system with pixel sizes of 172 for MIPS-24, roughly half of the original pixel size 
of 275. Pixel sizes of the final mosaics at 70 and 160/im are also approximately half of the 
originals, i.e. 4'.'95 at 70/im and 870 at 160/im. Unlike the IRAC coadds, we incorporated the 
sigma- weighting algorithm because we found it gave less noisy MIPS mosaics than without. 
Standard linear interpolation was used in all cases. In all three channels, our post-pipeline 
MOPEX mosaics are superior to those provided by the automated post-pipeline reduction. 
The final coadded images are shown in Fig. 6 as contours superimposed on the IRAC 4.5 
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/im image. We have performed aperture photometry on the MIPS images in the same way 
as on the IRAC images. The growth curves for the three bands (filled circles for the 24 jum 
band, open squares for the 70 /iin band and crosses for 160 /im band) are given in the right 
panel of Fig. 4. The dashed lines represent the point-response function in each band. The 
24 /iin, 70 /im and 160 /im total fluxes are listed in Table 1. 

2.3. IRS Spectra 

Spectroscopy was performed in the staring mode with the Short Low module in both 
orders (SL1, SL2) and with both Short and Long High-resolution modules (SH, LH). The 
data consist of low-resolution spectra with a wavelength range from 5.2 to 14.5 /im and a 
spectral resolution of ~64-128, and of high-resolution spectra with a wavelength range from 
9.6 to 37.2 /im and a spectral resolution of ~600. Total integration times of 30sx4 cycles were 
obtained for SH, 14s x 8 for LH, 6s x 6 for SL1, and 6s x 6 for SL2. The central region of Mrk 
996 was centered in the slits by peaking up on a nearby 2MASS star. Individual bed frames 
have been processed by the S13.2.0 version of the SSC pipeline, which provides ramp fitting, 
dark current subtraction, droop and linearity corrections, flat-fielding, and wavelength and 
flux calibrations 2 . 

However, the pipeline does not include background subtraction. Hence, for the low- 
resolution spectra, we constructed a coadded background frame from the bed observations 
with the source in the opposing nod and off-order positions (see also Weedman et al. 2005). 
For the high-resolution SH and LH spectra, a 2D background image could not be constructed 
because of the small size of the slit. Therefore we subtracted the background from the SH 
and LH observations using the one-dimensional (ID) spectra. The level of the background 
was determined by minimizing the difference between the SL and SH+LH spectra over their 
overlap region (<1 /im), and the shape of the background was given by Spot. 

We extracted the source spectra with SPICE, the post-pipeline IRS package provided 
by the SSC. To maximize the calibration accuracy, the automatic point-source extraction 
window was used for all modules. For the SL spectra, this gives a 4-pixel (7'/2) length at 
6/im, and an 8-pixel one (14'.'4) at 12 /im; the slit width is 3'.'6 for both SL modules. At high 
resolution, the SPICE extraction is performed over the entire slit (4 / . / 7xll'. / 3 SH; ll"lx22"3 
LH). The individual spectra were then box-car smoothed to a resolution element, and clipped 
in order to eliminate any remaining spikes in the high-resolution data. The spectra were also 
examined for bad pixel removal. Finally, the two spectra for each module (one for each nod 
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position) were averaged. The final averaged IRS spectra obtained in the short-wavelength 
low-resolution mode (SL) and in both short and long wavelength high-resolution modes (SH, 
LH) are shown in Fig. 7. A blow-up of the spectral region around the [O iv] line is shown 
in Fig. 8. 

3. Complementary Observations 

3.1. UKIRT Imaging 

We have also acquired J (1.2 /im), H (1.6 /im), and K (2.2 /im) images of Mrk996, 
with the 3.8 m UKIRT1 equipped with IRCAM3 as part of our ongoing project of NIR 
imaging and spectroscopy of BCDs. The IRCAM3 plate scale is 0'.'28 pixel -1 , with a total 
field of view of 72'.'0 x 72 / . / 0. Source and empty sky positions were alternated, beginning and 
ending each observing sequence with a sky position. Before the beginning of each sequence, 
dark exposures were acquired with the same parameters as the subsequent science frames. 
Total on-source integration times were 1040s in K, 480s in H, and 240s in J. Individual 
frames were first dark-subtracted and then flat-fielded with the average of adjacent empty 
sky frames, after editing them for stars (to avoid "holes" in the reduced frames) and applying 
an average a clipping algorithm. The reduced frames were then aligned and averaged. All 
data reduction was carried out in the IRAF environment. 

Photometric calibration was performed by observing standard stars from the UKIRT 
Faint Standard List (Hawarden et al. 2001) before and after the source observations. Each 
standard star was measured in several different positions on the array and flat-fielded by 
dividing the clipped mean of the remaining frames in the sequence. To correct the standard- 
star photometry for atmospheric extinction, we used the UKIRT mean extinction coefficients 
of 0.102, 0.059, and 0.088 mag per air mass for J, H, and K, respectively. Formal photometric 
accuracy, as measured by the dispersion in the standard-star magnitudes, is 0.025 mag in J 
and H and 0.04 mag in K. Color images were derived by registering the images to a few 
tenths of a pixel with a cross-correlation algorithm, then subtracting the magnitude images. 

As for IRAC, we have fitted ellipses to the J, H and K images of Mrk996, and J, 
H, and K elliptically-averaged radial profiles have been derived. The K surface brightness 
profile, which reaches a limiting surface brightness of ~ 24 mag arcsec -2 , is shown in the 
top right panel of Fig. 5. The lower three panels show the J — H, H — K and J — K radial 
color profiles, constructed from the J, H and K surface brightness profiles. In the last three 
panels, the dotted line indicates the average color of the underlying disk. 
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3.2. MMT Optical Observations 

As the HST optical spectra (Thuan et al. 1996) have somewhat low spectral resolution 
(~ 6A in the blue or ~ 350 km s -1 ), higher resolution optical spectra are needed to allow 
for a deconvolution of the line profiles into narrow and broad components and better study 
the gas motions in the inner and outer regions. A new high signal-to-noise ratio optical 
spectrum of Mrk 996 was obtained with the 6.5 m MMT on the night of 2006 December 
15. Observations were made with the Blue Channel of the MMT spectrograph. We used a 
l'.'5xl80" slit and a 800 grooves/mm grating in first order. The above instrumental set-up 
gave a spatial scale along the slit of 0'/6 pixel -1 , a scale perpendicular to the slit of 0.75A 
pixel -1 , a spectral range 3200 - 5000A and a spectral resolution of 3A (FWHM). The slit 
was oriented along the parallactic angle and the total exposure time was 25 minutes. The 
observations were broken up into 3 subexposures to avoid saturation of the brightest lines. 
The seeing was about 0'/8. The Kitt Peak IRS spectroscopic standard stars Feige 110 and 
G191B2B were observed for flux calibration. Spectra of He-Ar comparison arcs were obtained 
after the observations to calibrate the wavelength scale. 

The data reduction procedures are the same as described in Thuan & Izotov (2005). The 
two-dimensional spectra were bias subtracted and flat-field corrected using IRAF. We then 
use the IRAF software routines IDENTIFY, REIDENTIFY, FITCOORD, TRANSFORM 
to perform wavelength calibration and correct for distortion and tilt for each frame. Night 
sky subtraction was performed using the routine BACKGROUND. The level of night sky 
emission was determined from the closest regions to the galaxy that are free of galaxian 
stellar and nebular line emission, as well as of emission from foreground and background 
sources. A one-dimensional spectrum was then extracted from the two-dimensional frame 
using the APALL routine. Before extraction, the three distinct two-dimensional spectra 
of Mrk 996 were carefully aligned using the spatial locations of the brightest part in each 
spectrum, so that spectra were extracted at the same positions in all subexposures. We then 
summed the individual spectra from each subexposure after manual removal of the cosmic 
rays hits. The spectra obtained from each subexposure were also checked for cosmic rays 
hits at the location of strong emission lines, but none were found. 

Particular attention was paid to the derivation of the sensitivity curve. It was obtained 
by fitting with a high-order polynomial the observed spectral energy distribution of the two 
standard stars. Because the spectra of these stars have only a small number of a relatively 
weak absorption features, their spectral energy distributions are known with good accuracy 
(< 1%). Moreover, the response function of the CCD detector is smooth, so we could derive 
a sensitivity curve with an accuracy better than 1% over the whole blue optical range. 

The resulting spectrum of Mrk 996 is shown in Fig. 9. The spectrum is very rich in 
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emission lines, which we have identified and labeled. Because of its higher spectral resolution 
and signal-to-noise ratio, the MMT spectrum shows many more weak emission lines than the 
HST spectrum, confirming, strenthening and extending many of the findings of Thuan et al. 
(1996): 1) the He I line intensities such as that of the He I A4471 emission line are unusually 
large. There are many weak He I lines that are not ordinarily seen such as the A3820, 
A4121, A4143 and A4388 lines; 2) the widths of the emission lines of high-ionization stages 
are broader than those of low-ionization stages; 3) there are strong, broad Wolf-Rayet bumps 
at A4640 and A4686, suggesting the presence of a large WN stellar population. Thuan et al. 
(1996) have also detected a C IV WR bump at A5808, indicative of WC stars. In addition, 
the spectrum shows the presence of numerous weak [Fe ii] and [Fe ill] lines. The spectrum 
goes far into the blue, which allows us to check for the presence of the [Ne v] A3426 line, an 
indicator of hard ionizing radiation since its ionization potential is 7.1 Rydberg. It is not 
seen. The [Fe v] A4227 with an ionization potential of 4 Rydberg is also absent. The He II 
A4686, also with an ionization potential of 4 Rydberg, is seen but its broad width suggests 
that it originates in stellar winds of WR stars, not from the ionized interstellar gas. The 
Balmer jump, clearly seen at A3646, is useful for determining the temperature of the H + 
zone (Guseva, Izotov & Thuan 2006). 

The observed line fluxes F(X) normalized to F(H/3) and multiplied by 100 and their 
errors, are given in Table 2. They were measured using the IRAF SPLOT routine. The 
line flux errors listed include statistical errors derived with SPLOT from non flux-calibrated 
spectra, in addition to errors introduced in the standard star absolute flux calibration. Since 
the differences between the response curves derived for the two standard stars are not greater 
than 1%, we set the errors in flux calibration to 1% of the line fluxes. The line fluxes were 
corrected for both reddening (Whitford 1958) and underlying hydrogen stellar absorption 
derived simultaneously by an iterative procedure as described in Izotov, Thuan & Lipovetsky 
(1994). The corrected line fluxes 100 x I (X)/ 1(11/3) and equivalent widths EW(A) are also 
given in Table 2. 

Comparison of the relative line intensities given in Table 2 with those measured by 
Thuan et al. (1996) shows general good agreement. The largest deviation concerns the [O 
n] A3727 line which is about 1.5 times larger. We attribute this difference to an aperture 
effect: the HST aperture is circular with a diameter of / /86 while the MMT slit width is l'/5. 
Since the low-ionization ions are produced in the outer less dense regions, a larger aperture 
would have a larger [O n] to B./3 ratio. 
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4. Imaging Results: The Infrared Morphology of Mrk 996 

We now use the UKIRT near-infrared and the Spitzer IRAC and MIPS mid-infrared 
imaging data to study the origin of the infrared emission in Mrk 996. 

4.1. An exponential disk of > 1 Gyr old stars 

The top left and right panels of Fig. 5 show respectively the IRAC and fT-band surface 
brightness profiles of Mrk 996, which are very similar to the V and I profiles given by Thuan 
et al. (1996). In the radius range r > 4", they can be fit by an exponential law of the form 
Hk = 17.1 + 0.29 r("), /X4.5 = 16.7 + 0.25 r("), corresponding to scale lengths a K x ~ 
~0.4kpc, similar to the V and / scale lengths. This means that for r > 4", there is no 
V — K color gradient, and hence no stellar population change in the underlying extended 
low surface brightness component. This is also consistent with the constant V — I (= 0.80) 
(Thuan et al. 1996), J -H (=0.57) and H - K (=0.14) colors observed for r > 4" (panels 
2 and 3 of Fig. 5). These colors have been corrected for Galactic extinction Ay=0.146mag 
(Schlegel et al. 1998) with the Cardelli et al. (1989) interstellar extinction curve. 

The colors correspond to a coeval stellar population formed about 1 Gyr ago according 
to the models of Vazquez & Leitherer (2005). This age is fairly well constrained by the 
models for the optical colors, but in the NIR, models are still unable to completely predict 
broadband colors for intermediate-age stellar populations at sub-solar metallicities (e.g., 
Origlia et al. 1999; Mouhcine & Langon 2002; Vazquez & Leitherer 2005). The influence of 
red supergiants and asymptotic giant branch stars at ages of <1 Gyr increase with decreasing 
metallicities, and predicted NIR colors tend to be too blue in J — H and too red in H — K 
as compared to the observed ones. In any case, the observed NIR colors are clearly not 
those of old stellar populations in spiral disks, dominated by low- mass giants (e.g., de Jong 
1996; Peletier & Balcells 1997; Hunt et al. 1997). The oldest stars are probably older than 
1 Gyr as the star formation in Mrk 996 is likely more spread out in time as compared to the 
instantaneous star formation models we have used. 



4.2. A very red nucleus 

Just like the V and / surface brightness profiles (Thuan et al. 1996), the K profile shows 
a very steep gradient toward the nucleus. While the V — I profile becomes rather blue, in 
the near-infrared the striking feature of the nuclear region is its extremely red J — K color, 
equal to £ 1.8, as compared to the J — K color of the disk equal only to ~0.7. Examination 
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of the lower right panels of Fig. 5 shows that the redness of the nuclear region comes mainly 
from its H — K color (equal to ~1.0 as compared to the H — K = 0.1 of the disk), and 
less from its J — H color. The red colors cannot be produced by gaseous emission as pure 
gaseous colors are J — H = 0.0 and H — K = 0.6 (Hunt, Thuan, & Izotov 2003). The redness 
is more important at longer wavelengths, suggesting that it may be caused by a very hot 
dust component in the nucleus of Mrk996, similar to, although more extreme than, other 
star- forming galaxies (Hunt et al. 2002). 

Since the brightness gradient toward the nucleus is very steep, we have checked whether 
variable seeing for different exposures can be responsible for the very red colors observed. We 
found that the red J — K color of >2 shown in the color profile is slightly enhanced because 
of the better seeing in the K band image than in the J band image. However, the very red 
H — K color cannot be attributed to seeing effects. The seeing was 0'.'9 (FWHM) during the 
H band observations, and 0'/8 in the K band observations. We have convolved the -fT-band 
image to the if-band seeing, and produced a new H — K color image. The nuclear color 
in this new image is ~1.0, just as in the profile shown in Fig. 5. We have also checked the 
J — H color. The registration process necessary to align the images already degrades the H 
band image to the seeing of the J-band image of l'/O. The resulting nuclear J — H color is 
~0.8. Thus, the nuclear colors are still very red after correction for seeing effects: H — K = 
1.0, and J-H = 0.8. 

Fig. 10 shows in a ( J — H , H — K) diagram, as filled rectangles connected by a dotted 
line and from right to left, the colors of the nucleus, of a circumnuclear location, and of a 
disk location in Mrk996, together with those of IR-bright starburst galaxies (Hunt et al. 
2002) (open stars), of selected BCDs (asterisks) including NGC 5253, SBS 0335-052 (Hunt 
et al. 2002), and Haro 3 (Johnson et al. 2004), and of Seyfert galaxies (filled circles) (Alonso- 
Herrero et al. 2001). We have also shown in Fig. 10 mixing curves that illustrate how adding 
various physical components to the stellar component modify the JHK colors of the disk of 
Mrk996, corrected for the K effect and Galactic extinction. The five solid lines from top to 
bottom show respectively the effects of dust extinction, of a hot dust component at 600 K and 
1000 K, of free-free gas emission and of an A-star population. It can be seen that the JHK 
colors of the disk of Mrk996 are similar to the colors of other BCDs, and slightly bluer than 
those of normal (metal-rich) starbursts, but that the H — K color of the nucleus of Mrk996 
is considerably redder. While other BCDs show a clear contribution from ionized gas, the 
JHK colors of Mrk996 show a significant component of hot dust, with a dust temperature 
between 600 and 1000 K. 

What is the appearance of the nucleus at longer wavelengths? Examination of the MIPS 
images (Fig. 6) shows that it is clearly detected at 24, 70 and 160 /iin. The emission at 24 
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and 70 /im is virtually point-like, as shown by the characteristic Airy rings seen in the top 
two panels of Fig. 6, and from the growth curves in Fig. 4. However, the emission at 
160 /im appears to be more extended, as is also evident from its growth curve, suggesting a 
significant contribution from the galaxy disk. 

4.3. The circumnuclear PAH emission 

Examination of Fig. 1 and of the growth curves shows that both the 4.5 and 8.0 jum 
emissions are more extended than that of a point source. While the 8 /im emission comes 
mainly from the nuclear and inner-disk regions (out to a radius of ~10" or 1.1 kpc), the 
4.5 /xm image shows extended emission out to a radius of ~ 20" or 2.2 kpc, comparable in 
extent to the optical light (see the left panel of Fig.2). Because of the relative lack of 8 /im 
emission in the outer disk, the flux ratio map (Fig. 3) is very noisy there. Examination of 
the right panel Fig.2 shows that the 4.5/iin peak coincides with the bright blue star-forming 
nucleus of Mrk996 rather than with the red dust patch to the northwest (see Fig. 4a of 
Thuan et al. 1996). This suggests that the MIR emission is mainly associated with the 
source responsible for the very red nuclear J — K color. Fig. 3 shows that the 4.5/8 flux 
ratio is relatively constant (~0.7) outside the nuclear region (see §2.1), although beyond r 
~18", the color map becomes too noisy to give reliable colors. This can also be seen in the 
lower left panel of Fig. 5 where, for radii larger than ~ 6", the radial color profile shows a 
relatively constant [4.5] — [8.0] color of 1.5, corresponding to a 4.5/8 flux ratio of ~0.7. The 
apparent gradient in the circumnuclear region at the radius of ~ 4" is due to the mismatch 
of the diffraction-limited IRAC point-response functions at 4.5 /im and 8 /im. Because the 
4.5/8 color of stellar populations is ~ 2.5, roughly independently of age and metallicity 
(Leitherer et al. 1999), it is clear that the infrared emission of Mrk996 is dominated by the 
ISM, and not by stars. 

If we interpret the 8 /im emission in the usual way and attribute it primarily to PAHs, 
then the 4.5/8 ratio ~0.7 implies a contribution from PAHs in the ISM around the nuclear 
star-forming region of Mrk996. This interpretation is indeed supported by the detection of 
strong PAH features in the IRS spectrum discussed in the next section. These PAHs are 
probably located in the outer low-density region with 400 pc < r < 1000 pc, since they 
are likely destroyed by the strong UV radiation field in the nuclear region. However, our 
observations do not have sufficient spatial resolution to verify this. 
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5. Spectroscopic results 

5.1. PAH features 

The IRS spectra (Fig. 7) show clearly PAH features at 5.7, 6.2, 7.7, 8.6, 11.2 and 12.8 
/im. The flux, the equivalent width EW and the full width at half maximum FWHM of each 
PAH feature have been derived by fitting Lorentzian profiles (Gaussian profiles do not give 
good fits) with the SPLOT task in IRAF. They are given in Table 3. The continuum was 
linearly interpolated in two sections, one from 6-9 /im, and the other from 10.5-13.9 /iin. Stan- 
dard deviations of repeated measurements are given in parentheses. The relative strengths 
of the three main PAH emission features can be used as diagnostics to identify PAH sizes 
and infer the neutral-to- ionized gas ratio (Draine & Li 2001). According to the models of 
Draine & Li (2001), the observed PAH flux ratios f n . 2//V.7 equal to 0.40 and f^l h.i equal 
to 0.27 indicate that the PAHs in Mrk996 are predominantly neutral and small, containing 
a few hundred carbon atoms, similar to the normal galaxies studied by Helou et al. (2000). 
In other words, just as Haro3, Mrk996 lies in the "normal galaxy" region in the PAH di- 
agnostic diagram rather than in the starburst galaxy region which includes objects such as 
M82 and NGC253, that are dominated by ionized PAHs. This would imply, at face value, 
that the PAHs in Mrk996 arise mainly from the normal cold and neutral ISM surrounding 
the nuclear region, but not from the star-forming region itself. However, the PAH EWs 
are generally high, ~ 3.6/im for the 7.7 /im line, which is more typical of starburst galaxies 
(Brandl et al. 2004) than of BCDs. Wu et al. (2006) found that the EW of the 7.7 /im line 
of all the BCDs in their sample to be less than 1 /im, considerably smaller than the value 
of 3.6 /mi in Mrk996. For the metallicity of Mrk996, the EWs of its 6.2 and 11.2 /im PAH 
features exceed by more one order of magnitude the values given by the mean EW (PAH) - 
metallicity relations given by Wu et al. (2006) for their BCD sample. 

Just as for Haro3, we find that the PAH emission features in Mrk996 are relatively 
narrow. They are narrower than those predicted by the model for the diffuse Galactic ISM 
at high latitudes (Li & Draine 2001). While the 6.2 /im and 7.7 /im features have FWHMs 
similar to those of the starburst galaxies M 82 and NGC 253, the 8.6 /im feature is about half 
as wide while the 11.3 /tm feature is twice as wide. The best matches to the FWHMs of the 
PAH features in Mrk996 appear to be those of the reflection nebulae region NGC 2023. Thus, 
neither the diffuse Galactic ISM model nor the starburst model of Li & Draine (2001) and 
Draine & Li (2001) can account fully for all the observed properties of the PAHs in Mrk996. 
Even though the IRS spectrum of Mrk996 may contain a conspicuous ISM component, 
Galactic ISM spectra are not necessarily good templates for all PAH emission in this object. 
This suggests that there are many factors that determine the properties of PAH features 
(e.g., chemical abundance, ISM energetics, etc.) which are not yet taken into account fully 
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in the models. A more definitive observational test of the models, with a fuller exploration 
of the parameter space, will be possible with our larger sample of 23 BCDs. 

5.2. Infrared Fine-Structure Lines 

The spectrum in Fig. 7 also shows several fine-structure lines. The [Siv]A10.51, 
[Nell]A12.81, [Nem]A15.55, [Sm]A18.71, 33.58, and [Oiv]A25.89 fim are seen. In the same 
way as for the PAH features, we used SPLOT to fit the IR fine-structure emission lines and 
obtain fluxes and EWs. A deblending procedure was adopted to accurately measure the 
emission lines at wavelengths near PAH features (e.g., [Nell]). The lines were best fit with 
Gaussian profiles. A single continuum was linearly interpolated from 10.5 to ~ 15 /im. For 
the longer wavelengths, the local continuum was fit by linearly interpolating over adjacent 
line-free regions. The continuum la uncertainty in the SH module is ~ 0.18 x 10~ 17 Wm~ 2 , 
and ~ 0.11 x 10~ 17 Wm~ 2 in LH. The fluxes and other parameters resulting from the fits 
are given in Table 4. 

Several fairly high-ionization lines are detected in Mrk996. The [O lv]A25.89/xm line 
with an ionization potential of 54.9 eV, just beyond the He II edge at 54.4 eV, is faint, but 
present at the 4<r level. A blown-up view of the spectral region around the [O iv] line is 
shown in Fig. 8. However its presence raises a puzzle. As discussed before, the optical 
spectrum of Mrk 996 does not show the presence of the high-ionization [Fe v] A4227 emission 
line which also has an ionization potential of 54.4 eV. The He II A4686 emission line, also 
with an ionization potential of 54.4 eV, is seen, but its broad width suggests that it originates 
in WR stellar winds rather than in the ionized interstellar gas. This implies that the ionizing 
radiation in the optical line-emitting region is less hard than that in the MIR line-emitting 
region. 

We can have some idea of the hardness of the radiation in the MIR-emitting region by 
using line ratios of different ionic species of the same element, such as [Nelll] and [Nell], 
since these are sensitive to the shape of the spectrum of the ionizing radiation field. The ratio 
F([Ne ill] )/F( [Nell]) is 2.8, placing Mrk 996 near the high-excitation end of the starbursts 
studied by Verma et al. (2003), similarly to NGC5253 and Haro3. Because PAHs tend to 
be depleted in hard radiation fields, their EW would be expected to be negatively correlated 
with the [Ne ill] /[Nell] line ratio. Such a trend was indeed found by Wu et al. (2006) for 
their BCD sample. The anticorrelation between the PAH strength and the hardness and 
luminosity of the UV radiation field was also discussed by Beirao et al. (2006) in the case 
the irregular galaxy NGC 5253 and by Lebouteiller et al. (2007) in the case of the massive 
young cluster NGC 3603. However, as discussed before, compared to the Wu et al. (2006) 



- 15 - 



correlation, the PAHs in Mrk996 have considerably higher EWs for their metallicity. This 
may again be due to a significant PAH component from the general ISM in Mrk996 as 
compared to the BCDs studied by Wu et al. (2006), as the region of Mrk996 covered by the 
IRS slit is relatively large, about 700 pc. The large EWs may also be the results of different 
measurement techniques, as Wu et al. (2006) fit a spline to the local underlying continuum, 
while we linearly interpolate over two wide sections of continua. 

The long-wavelength IRS spectrum (Fig. 7) shows an emission feature at A29.84 fim, 
the wavelength of one of the strongest H 2 lines in this spectral region. Such a detection 
would not perhaps be surprising, as water can form in shock-heated gas associated with 
outflows (e.g., Melnick et al. 2008, and references therein). However, we would also expect 
to detect a wealth of H 2 lines which are not clearly seen in our spectrum (Fig. 7). Hence, it 
is not clear whether this is a real feature or not. We will investigate the frequency of water 
emission in a future paper on our larger BCD sample. 

5.3. The Infrared Spectral Energy Distribution 

We combine the available photometric and spectroscopic data to construct the IR spec- 
tral energy distribution (SED) in Fig. 11. We have used the total emission from Mrk996, 
integrated over the entire galaxy including the disk. The 8 /iin IRAC point which is higher 
than the corresponding IRS point suggests that a significant part of the PAH emission in the 
disk lies in a region more extended than the one covered by the IRS SL slit. The coincidence 
of the MIPS 24 /xm point with the IRS spectral one lends confidence to our data reduction 
procedures, including background subtraction and flux calibration. 

It can be seen that the SED has a broad peak from 80 to 160 fxm, implying that the dust 
in the disk of Mrk996 must be rather cool, <20K. Indeed, this broad peak toward longer 
wavelengths is almost certainly due to the "mixed morphology" of Mrk996. Up to <70 fim, 
Mrk996 is virtually a point source, dominated by the compact nuclear star-forming region; 
but at longer wavelengths, the disk dominates with its contribution of cooler dust. Hence, 
the SED does not show a well-defined spectral peak, but rather a smeared-out one, due to 
the contributions of dust at different temperatures from different regions. 

We have compared in Fig. 11 the SED of Mrk996 with model SEDs of starburst 
galaxies such as Arp220 and M82, taken from Silva et al. (1998). These models, normalized 
to continuum emission between 20 and 30 /im, clearly do not fit the SED of Mrk996 (see also 
the discussion of Hunt et al. 2005): Mrk996 shows no silicate absorption feature at ~10/xm, 
relatively little PAH emission compared to a luminous solar-metallicity starburst galaxy, 
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and a very different MIR continuum slope. Moreover, in Mrk996, the unabsorbed emission 
that emerges at NIR and optical wavelengths is greater relative to the dust emission, when 
compared to the prototypical starburst galaxies M82 and Arp220. This is most probably 
due to a lower extinction and thus to a lower fraction of dust reprocessing, but could also 
be due to differing dust content relative to the stars. 

We have also compared the SED of Mrk996 with that of the Seyfert 1 galaxy NGC4151 
(Buchanan et al. 2006), normalized over the 9 -15 jum region (Fig. 12). While this SED 
gives low residuals in the normalization zone, il fails to fit the SED for wavelengths longer 
than 15 /im. Clearly, dust reprocessing of the H II region radiation dominates the SED of 
Mrk996 at longer wavelengths. The surrounding cooler ISM also contributes significantly, 
resulting in a roughly constant (in flux density) SED from 70 to 160 /im. 

In summary, like most other star-forming galaxies (see the discussion of the BCD Haro3 
in Hunt et al. 2006), Mrk996 is a composite entity in the IR; different components and mor- 
phologies are seen at different wavelengths. In Mrk996, at JHK we see the clear signature 
of a very red nucleus, and the evolved extended underlying stellar population in the disk; 
similarly, at 4.5 /im we see a nuclear hot dust continuum and the evolved stars in the disk; 
at 8/iin, we see mostly extended PAH emission in the general ISM; at 24 and 70 /im, we 
see unresolved emission from warm dust associated with the nucleus; and at 160 cooler 
extended dust emission in the disk. 



6. CLOUDY modeling of the emission-line fluxes 

6.1. The narrow and broad line components of Mrk996 

Thanks to a higher spectral resolution than that of the IRS spectrum, it can be seen 
clearly from the MMT optical spectrum (Fig. 9) that some lines have two components: a 
narrow component at high intensity levels and a broad component at low intensity levels. We 
have used SPLOT to deconvolve the profiles of these lines into narrow and broad components. 
The deconvolution results are shown in Table 5. It gives for both components the line flux of 
each detected line normalized to the H/3 flux and its FWHM. The narrow component of the 
[O ill] A4363 line is negligible. Three facts are to be noticed in Table 5. First, the narrow 
line component has FWHMs varying between 2.6 and 3.5 A, comparable to the spectral 
resolution. It is thus not resolved, giving upper limits of gas velocities between 87 and 135 
km s -1 . On the other hand, the broad line component has FWHMs varying between 6.1 and 
8.5 A, corresponding to gas velocities between 235 and 330 km s" 1 . Second, the ratio of the 
auroral [O ill] A4363 line to that of the nebular [O ill] A4959 line for the broad component, 
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equal to ~ 0.45, is unusually high compared to its usual value of ~ 0.1 for an electron 
temperature of 20 000 K, suggesting collisional deexcitation of the [O ill] A4959 line. Third, 
the intensities for the He I lines relative to B./3 are different for the narrow and broad line 
components. For the He I A3889 line, the relative intensity of the broad component is smaller 
than that of the narrow component by a factor of 1.4, while for the He I A4471 line, the 
reverse is true, the broad component being more intense by a factor of 2. 

Thuan et al. (1996) have suggested that these observational features can be understood 
if the optical lines arise in two main regions: (1) the broad lines in a very dense inner 
region (r < 100 pc) with large gas mass motions, most likely powered by stellar winds of the 
Wolf-Rayet stars, the presence of which is clearly indicated by the WR bump in the MMT 
spectrum (Fig. 9) , and (2) the narrow lines in a considerably less dense outer region (r > 
several hundred pc), with a density more comparable to those in normal H II regions, and 
smaller mass motions. In this type of model, the helium lines would arise in the dense inner 
region: the intensity of the He I A4471 emission line would be enhanced because of collisional 
excitation, while the intensity of the He I A3889 line which is sensitive to optical depth effects, 
would be decreased because of the higher optical depth. Thuan et al. (1996) were indeed 
able to construct an inhomogeneous two-zone CLOUDY model which accounts well for the 
observed line intensities of the UV and optical lines. Their best model has an inner region 
with a central electron density N e equal to 4.5 x 10 6 cm~ 3 and a density decreasing with 
radius r as r~ 2 out to a radius of 100 pc, and an outer region with a constant density N e = 
450 cm' 3 for radii between 100 pc and 320 pc. 

In addition to the presence of distinct broad and narrow line regions, there is also a trend 
of increasing line width with ionization potential, similar to Seyfert galaxies (De Robertis & 
Osterbrock 1984; Whittle 1985). However the gas clouds motions and densities in the broad- 
line region are not as large as in the case of an AGN. Also, the UV and optical spectra show 
no evidence of the very hard ionizing radiation that is usually present in an AGN. But, as 
discussed before, the detection of the [O iv] A25.9/xm emission line in the IRS spectrum (Fig. 
7) does imply the presence of hard EUV radiation in the MIR-emitting region. However, 
hard ionizing radiation, as seen below, can be produced by shocks or by WNE-w stars, not 
necessarily by an AGN. 

6.2. A hidden H II region 

We now use CLOUDY to check whether the two-zone model constructed by Thuan et 
al. (1996) to account for the UV and optical line fluxes observed by HST can also explain 
the MIR line fluxes observed by Spitzer. 
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The parameters of the CLOUDY model used here are very similar to those of model 
1 of Thuan et al. (1996) (the one with low nitrogen abundance). The slight differences are 
due to the updating of the CLOUDY code with new atomic parameters. However, since the 
Spitzer and MMT spectra have been obtained through significantly larger apertures than 
the HST aperture, we have increased the radius of the modeled H II region to ~ 580 pc in 
order to take into account the total Ha emission of Mrk996. This emission extends over 
an angular radius of ~ 5" (~520pc) (Gil de Paz et al. 2003), and is covered entirely by 
the LH IRS slits and mostly by the SH IRS slits. The model gives a H/3 luminosity L(R(3) 
= 2.34xl0 40 erg s _1 , corresponding to a number of ionizing photons per second, log Q(R) 
= 52.75. At the distance of 21.7 Mpc, this luminosity corresponds to a H/5 flux /(H/3) = 
4.18xl0~ 16 W m" 2 . This flux is very similar to the B./3 flux of 4.40 xl0~ 16 W m~ 2 derived 
from the extinction-corrected Ha flux measured for Mrk 996 by Gil de Paz et al. (2003), 
adopting an extinction coefficient C(H/3) = 0.53 and a Ha/H/3 ratio of 2.8. As for the run 
of the number density with radius r, we have also adopted for the inner part (1 pc < r < 
100 pc) a r~ 2 law, with a maximum log N e = 6.66 at the inner boundary (r = 1 pc), and 
a constant log N e = 2.645 for r > 100 pc. The filling factor, log / = —4.58, characteristic 
of a very clumpy interstellar medium, is the same as in Thuan et al. (1996). Hereafter, this 
model will be referred to as Model I. 



6.2.1. Comparison of observed and modeled optical line fluxes 

First, we check that Model I is indeed able to reproduce the fluxes of the optical emission 
lines obtained from the MMT spectrum. Table 6 lists the extinction-corrected observed fluxes 
and Model I fluxes of several strong optical lines along with their ratios. It is seen that the 
observed high-ionization [O ill] and [Ne ill] emission lines are well reproduced by the model, 
but that the observed fluxes of of the lower ionization [O il] A3727 and H/3 emission lines are 
significantly smaller (by factors of 2 and 1.7 respectively) than the modeled ones. Evidently, 
the difference is due to aperture effects because the aperture of the MMT spectrum does not 
cover all the emission from the H II region, and the lower-ionization lines are produced further 
away from the nucleus than the high-ionization ones. Therefore, we have also calculated 
Model la, which has the same parameters as Model I, except that the radius of the modeled 
H II region is equal to 410 pc instead of 580 pc. This radius was chosen so that the modeled 
H/3 flux matches the observed H/5 flux within the MMT aperture. Compared to Model I 
fluxes, the fluxes of the high-ionization lines in Model la remain unchanged, while the flux 
of the [O il] A3727 emission line is considerably lower, bringing it into good agreement with 
the observed one, as shown in the last column of Table 6. Thus, we conclude that CLOUDY 
models with the parameters adopted above do indeed reproduce well the observed fluxes of 
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the optical lines. 

6.2.2. Comparison of observed and modeled MIR line fluxes 

We next compare the observed MIR emission line fluxes with those predicted by Model 
I. Examination of columns 2 and 3 of Table 7 shows that, despite similar observed and 
modeled fluxes of the optical lines, Model I predicts significantly lower MIR line fluxes, as 
compared to the observed ones, except for [Ne n] A12.8 /im. There could be at least two 
explanations of this discrepancy: 1) aperture effects and 2) the presence of an additional 
obscured H II region that is not seen in the optical range, but that contributes to the MIR 
emission. 

We consider first the possible role of aperture effects. The angular radius of the optical 
H II region of ~ 5" is significantly greater than half of the slit width used in the optical 
observations. It is greater than or comparable to half of the slit width used for the IRS 
spectra. Therefore, aperture effects may be important. In particular, these effects may play 
a role for the hydrogen lines and the low-ionization forbidden optical [O n] A3727 and the 
MIR [Ne n] A12.8/zm lines because a significant fraction of emission in these lines comes from 
the outer zones of the H II region. For higher-ionization lines, the aperture effects should be 
smaller because those lines arise in the inner zones of the H II region. To check how aperture 
effects may influence the emission line fluxes, we show in Figure 13 the fluxes F predicted by 
CLOUDY of different optical and MIR lines emitted by a H II sphere with a varying radius r, 
keeping all other parameters equal to those of Model I. All fluxes are normalized to the fluxes 
F max from the sphere with the 580 pc outer radius of the H II region in Model I. We have also 
shown by vertical lines the "equivalent" radii of the different apertures used in the optical 
and MIR observations, defined as a/ S/tt, where S is the area of the aperture. We also mark 
the location of the outer radius of the H II region in Model I. The F/F max ratios cannot be 
used directly for aperture corrections since they are derived for spheres, while the observed 
emission comes from a cylindric column with equivalent radius r. They provide however 
upper limits to the aperture corrections and Fig. 13 can be used for a qualitative analysis. 
It is seen from that figure that aperture corrections are important for some optical lines, 
especially the [O il] A3727 line, as confirmed by comparison of the HST and MMT fluxes: 
the MMT flux (Table 2) is 1.5 times larger than the ffSTflux (Thuan et al. 1996). The HST 
and MMT B.f3 fluxes are respectively ~ 3 and 2 times lower than the one inferred from the 
Ha flux by Gil de Paz et al. (2003). On the other hand, the aperture correction is negligible 
for higher-ionization optical lines, such as [O ill] A5007. The aperture corrections for the 
MIR lines are significantly lower because of the larger apertures used. Fig. 13 shows that 
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the IRS SH aperture includes nearly all the H II region emission and the IRS LH aperture 
includes all of it. Thus, aperture effects cannot account for the large differences between 
the observed and modeled fluxes of the MIR lines, except perhaps for the [Ne n]A12.8/xm 
emission line which has the largest aperture correction. We conclude that, to account for 
the MIR emission, we need to postulate the existence of a hidden H II region in Mrk 996 
that contributes to the MIR emission, but not to the optical emission. 

We have found that the best CLOUDY model of this obscured H II region (hereafter 
Model II) is characterized by the following parameters: a number of ionizing photons log 
Q(H) = 52.425 - ~ 2 times lower than for Model I -, an effective temperature T e // — 50000 
K, - slightly higher than in Model I, suggesting a younger H II region -, a constant number 
density of 300 cm -3 , a filling factor log / = -2, an H II region radius of ~ 100 pc and the 
same chemical composition as in Model I. The predicted fluxes of the MIR lines for Model 
II are shown in Table 7. It is seen that, except for the [Ne n]A12.8/im line, the predicted 
fluxes of the MIR lines in Model II (the hidden H II region) are all significantly higher than 
those in Model I (the visible H II region), despite a number of ionizing photons which is 
~2 times lower. The main reason is that, in Model I, the infrared fine-structure line fluxes 
are strongly reduced by collisional deexcitation of the upper levels in the very dense central 
part of the visible H II region. At the critical densities, the line ratios will change as some 
lines will be suppressed less strongly as compared to other lines that will be suppressed more 
strongly. The critical number densities for collisional deexcitation of the detected MIR lines 
are in the range ~I0 3 - fO 5 cm -3 (Martin-Hernandez et al. 2002), lower than the central 
number density. Collisional deexcitation is however not important for the [Ne ll]A12.8/im 
line because it has a very high critical density of 6.1 xfO 5 cm~ 3 and in addition it arises in 
the outer less dense zones of the H II region in Model I. Adding the line fluxes from Models 
I and II (columns 3 and 4 of Table 7), we see that the modeled and observed fluxes of the 
MIR lines are in general good agreement, their ratios (column 8 of Table 7) being close to 
unity. How much dust extinction would it take to hide the H II region of Model II ? If we 
set the upper limit of its H/3 flux to be f 0% that of the visible region, then its Ay would be 
> 4 magnitudes. 



6.3. A hidden source of hard ionizing radiation 

There is, however, one glaring discrepancy in column 8 of Table 7: the model does not 
predict any [O iv] A25.9/zm emission, while this line is clearly detected in the IRS spectrum. 
The ionization potential of [O iv] is 54.9 eV, just beyond the He II edge at 54.4 eV. As 
discussed before, we have searched the optical spectrum of Mrk 996 for the high-ionization 
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[Fe v] A4227 and He II A4686 nebular emission lines which both have an ionization potential 
of 54.4 eV, but none were seen. This implies that the hard ionizing radiation comes entirely 
from a region that is optically hidden, which is not necessarily spatially coincident with the 
optically hidden H II region discussed before. 

Mrk 996 is not the only BCD to reveal harder ionizing radiation in the MIR than the op- 
tical. The same phenomenon occurs in the BCD Haro 3, the most metal-rich [12+log(0/H) = 
8.32] BCD of our Spitzer sample. Haro 3 also shows a well detected [O iv] A25.9/im emission 
line in its IRS spectrum, while the [Fe v] A4227 and He II A4686 nebular emission lines are 
conspicuously absent from its optical spectrum (Hunt et al. 2006). We will have to wait for 
the analysis of our whole Spitzer BCD sample to tell whether that is a general property of 
BCDs, or whether this occurs only in exceptional cases. 

What is the origin of the hard ionizing radiation beyond 4 Ryd in Mrk 996? Can stellar 
radiation be responsible? We have run Costar models by Schaerer & de Koter (1997) 
with stars having the highest effective temperature (T e ff = 53 000 K) and the hardest 
radiation possible (corresponding to the hottest 03 stars). We failed to reproduce the 
observed intensity of the [O iv] A25.9yum by a factor of several. We have then considered 
Wolf-Rayet stars of type WNE-w. According to calculations by Crowther et al. (1999), 
models for WNE-w stars (early nitrogen Wolf-Rayet stars with weak lines) show a strong 
ionizing flux above 4 Rydberg (A < 228 A), in contrast to the WCE and WNL stars which 
show negligible fluxes above 3-4 Rydberg. These stars have transparent winds, and therefore 
their Q(> 4 Kyd)/Q (> 1 Ryd), where Q is the ionizing flux, is some 100 times higher 
than the same ratio predicted for the hottest O stars by the Costar models (Schaerer & 
de Koter 1997). We have run CLOUDY models with the ionization spectrum of a WNE- 
w star (Crowther et al. 1999). We found that with a log Q(H) = 51.425, some 25 times 
less than that from the bright H II region, and some 10 times less than the Q(H) of the 
hidden H II region, the observed intensity of the [O iv] line can be reproduced (Table 7). 
The fraction of WNE-w to O stars is several percent if the same Q(H) (log Q(H) ~ 49) 
is assumed for single O and WNE-w stars. This estimate may be wrong by a factor of 10 
because of the uncertainty in the Q(H) of WNE-w stars. The WNE-w stellar spectrum does 
not produce much of the lower-ionization species, so that the agreement obtained previously 
between observations and models for the line intensities of the lower-ionization species is not 
destroyed. A weak [Ne v] A24.3yum line is predicted by CLOUDY, with an intensity some 9 
times weaker than that of the [O iv] line intensity, consistent with its non-detection. How 
plausible is the presence of such a WNE-w stellar population in the central star-forming 
region of Mrk 996? A ratio N(WNE)/N(0) of a few percent appears high because, at the 
metallicity of Mrk 996 it is roughly equal to the number ratio of WR stars of all types to 
O stars (Guseva, Izotov & Thuan 2000), and because WNE-w stars are rare among WR 
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stars. Moreover, from the equivalent width of H/3 (106.9 A), we estimate the age of the 
bright visible H II region to be 3-4 Myr. WNE stars are thought to result from mass transfer 
in close binary systems, so they dominate the stellar population at ages > 10 Myr, when 
O stars and WR stars of other types are gone. Thus the visible H II region is too young 
to contain WNE stars. If one wishes to invoke WNE stars as the source of hard radiation, 
one would have to postulate that these stars are not visible in the optical range. All things 
considered, we do not consider the WNE-w hypothesis very likely. 

We have also considered the case of a non-stellar source of ionizing radiation, that 
emitted by an accretion disk around an intermediate-mass black hole in the center of Mrk 996. 
Thus, we have also run CLOUDY models with an ionization spectrum in the form of a power- 
law: f(v) oc u^ 1 . We found that, in order to reproduce the line intensity of the [O iv] line, 
the number of ionizing photons coming from the AGN has to be log Q (H) = 51.125, some 50 
times smaller than that coming from the bright H II region. Again, the previous agreement 
between models and observations for the lower-ionization species is preserved. The results 
are also given in Table 7. Again, a weak [Ne v] A24.3/im line is predicted, with an intensity 
some 3 times weaker than that of the [O iv] line intensity. This line can be even weaker if 
the power-law spectrum has a slope steeper than —1. In any case, such a low- intensity line 
would be hard to detect. Because we have to postulate that the AGN is completely invisible 
in the optical, we do not favor the AGN hypothesis as the most probable one. 

Finally, we discuss a third possible source of hard ionizing radiation. Izotov et al. 
(2001), Izotov et al. (2004) and Thuan & Izotov (2005) have suggested that fast radiative 
shocks moving through a dense interstellar medium with N e ~100 cm -3 can be possible 
sources of photons with energy >54 eV, and be responsible for the [Ne v] A 3426 (ionization 
potential of 7.1 ryd), [Fe v] A 4227 (ionization potential of 4 ryd), and He II A 4686 (ionization 
potential of 4 ryd) emission they observed in some BCDs. The ionizing spectrum of such 
fast shocks has been computed by Dopita & Sutherland (1996) for a gas of solar metallicity 
and various shock velocities. Recently, this work has been extended to environments with 
lower metallicities by M. G. Allen et al. 2004 (in preparation). Thus, we have also run 
CLOUDY models with an ionizing spectrum taken from the "Mapping" III Shock model 
library (http : / /cdswww.u — strasbg.fr/ alien /mappings _pagel.html), and characterized 
by a metallicity equal to that of the Small Magellanic Cloud (closest to the metallicity of 
Mrk 996) and a shock velocity equal to 250 km s _1 . This velocity corresponds to the FWHM 
of ~ 7 A observed for the broad component of the emission lines in Table 5. The CLOUDY 
results are given in Table 7 for the shock model with a number of ionizing photons log Q(H) 
= 51.425, the same as that for the model with the WNE-w stars. It can be seen that the 
predicted line intensities in the shock model are very similar to those in the WNE-w model, 
and thus can account equally well for the intensity of the [O iv] line. Again a weak [Ne v] 
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A24.3/im line is predicted which would be undetectable. We favor the shock hypothesis as 
the most plausible explanation for the [O iv] line, as Mrk 996 contains many WR stars with 
outflowing stellar winds and, likely, supernova remnants which will no doubt cause radiative 
shocks. 



7. Summary 

We have acquired Spitzer MIR, UKIRT NIR and MMT optical observations of the blue 
compact dwarf galaxy Mrk 996 to study its gas, dust and stellar content. This BCD, with a 
metallicity of about 1/5 that of the Sun, has the peculiarity of possessing an extremely dense 
nuclear star- forming region: its central density is ~ 10 6 cm" 3 , some 4 orders of magnitude 
greater than the densities of normal H II regions. We have obtained the following results: 

(1) The nucleus of Mrk 996 is extremely red, with J — K = 1.8, and H — K = 1.0, probably 
due to very hot dust with a temperature between 600 and 1000 K. The optical spectrum of 
the BCD shows the high-ionization lines to have both broad and narrow line components, 
and a trend of increasing line width with increasing ionization potential. 

(2) The VIJHK colors of the underlying exponential disk are roughly consistent with the 
colors of a coeval stellar population with age > 1 Gyr. 

(3) Like most star-forming galaxies, Mrk 996 is a composite entity in the IR. We see extended 
photospheric emission from evolved stars, compact hot dust continuum coming from the 
nuclear star-forming region at 4.5 /xm, hot dust continuum and extended PAH emission 
coming mainly from the surrounding less dense ISM at 8 /im, compact small grain warm 
dust associated with the active star-forming nuclear region at 24 /im and 70 /im, and cooler 
extended dust emission associated with older stellar populations at 160 /im. 

(4) The IRS spectrum (Fig. 7) shows strong Polycyclic Aromatic Hydrocarbon (PAH) molec- 
ular emission, with features clearly detected at 5.7, 6.2, 7.7, 8.6, 11.2 and 12.7/im. The PAHs 
in Mrk 996 are predominantly neutral and small, similar to those found in normal spiral galax- 
ies, suggesting that they reside in the general ISM and not in the star-forming region. The 
PAH emission features are relatively narrow and their equivalent widths are generally high 
for the metallicity of Mrk 996, exceeding by more than one order of magnitude the values 
given for the mean EW(PAH)-metallicity relation derived by previous investigators. 

(5) Gaseous nebular line emission is seen. The IRS spectrum shows several fine-structure 
forbidden lines, including [Siv] A10.51, [Nell] A12.81, [Nelll] A15.55, [S in] A18.71, 33.48 
and [Oiv] A25.90/iin. 
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(6) We have used CLOUDY to model the line-emitting region. To account for both the 
optical and MIR lines, two H II regions are required: a) a very dense H II region that is 
seen in the optical range (Model I) and b) an optically obscured H II region (A v > 4 mag) 
with a constant number density of ~ 300 cm" 3 , typical of other H II regions (Model II). 
A two-zone model is required for the non-obscured H II region (Model I): a) a very dense 
nuclear region where the broad optical line components arise; from a central value of ~ 10 6 
cm" 3 , the density decreases with distance as r~ 2 until r ~ 100 pc; b) an outer zone for 100 
pc < r < 580 pc, with constant number density of ~ 400 cm" 3 . The density gradient is 
probably caused by large scale gas mass motions, powered by the stellar winds of Wolf-Rayet 
stars. 

(7) The UV and optical spectra show no evidence of the very hard ionizing radiation seen in 
AGN. The IRS spectrum does show however the presence of a faint [Oiv] line at 25.89 /im, 
indicating the presence of radiation as hard as 54.9 eV in Mrk996. This hard radiation is 
most likely due to fast radiative shocks in the ISM caused by outflowing stellar winds from 
WCE and WNL stars and/or by supernova remnants. A hidden AGN, or a population of 
hidden Wolf-Rayet stars of type WNE-w are are less likely sources of hard ionizing radiation. 

This work is based on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology, under NASA 
contract 1407. The MMT time was available thanks to a grant from the Frank Levinson Fund 
of the Peninsula Community Foundation to the Astronomy Department of the University of 
Virginia. We thank Jack Gallimore for providing us with the spectrum of NGC4151. We are 
grateful to the UKIRT Time Allocation Panel for a generous time allocation. Support for 
this work was provided by NASA through contract 1263707 issued by JPL/Caltech. T.X.T. 
and Y.I.I, also acknowledge partial financial support from NSF grant AST 02-05785. L.K.H. 
and Y.I.I, are grateful to the hospitality of the Astronomy Department of the University of 
Virginia where part of this work was done. 
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Fig. 1. — The 4.5/mi (IRAC2) image (left panel) and the 8.0/zm (IRAC4) image (right panel) 
of Mrk 996. The units given in the colorbar are MJy/sr, and the grey scale is logarithmic. The 
contours in both panels correspond to the 4.5/im image. The outermost contour corresponds 
to the sky level + 10a (a = 0.0028 MJy/sr). 
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Fig. 2.— The 4.5/im IRAC image of Mrk996 contoured on: (left panel) the HST/WFPC2 
F791W image (Thuan et al. 1996) and (right panel) the HST F569W-F791W color image 
(Thuan et al. 1996). The units in the color bar are counts (left) and magnitudes (right), and 
the grey scale is logarithmic. The contours in both panels correspond to the 4.5/im image 
(see Fig. 1). The outermost contour corresponds to the sky level + 10a (er = 0.0028 MJy/sr). 
The IRAC 4.5/xm brightness peak is centered on the blue nucleus (F569W— F791W ^ 0) of 
the galaxy, rather than on the red dust patch (F569W— F791W~3) to the northwest. 
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Fig. 3. — The 4.5/im IRAC image of Mrk996 contoured on the IRAC 4.5/8.0yum flux ratio 
image. The flux ratio units given in the color bar are in MJy/sr. As in Fig. 1, the outermost 
contour corresponds to the sky level + 10cr (er = 0.0028 MJy/sr). The local maximum of the 
4.5/8.0/xm flux ratio is on top of the 4.5/im brightness peak, which also coincides with the 
blue nucleus (F569W— F791W ^ 0) of the galaxy. There is very little 8 //m emission outside 
the circumnuclear region, so the color ratio map is considerably more noisy in the outer parts 
of the disk. 
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Fig. 4. — Growth curves for IRAC (left panel) and MIPS (right) photometry of Mrk996. In 
both panels, for each wavelength, the horizontal dotted line shows the adopted total flux of 
Mrk996, and the dashed line represents the point-source response function. In the left panel, 
the 4.5 /zm IRAC photometry is indicated by filled circles, and the 8.0 /zm photometry by 
filled squares. In the right panel, the MIPS-24/im photometry is shown by filled circles, the 
MIPS-70/im photometry by open squares, and the MIPS-160;um photometry by x. Mrk996 
is extended at shorter wavelengths and possibly at 160/im, but appears essentially point-like 
at 24 and 70/iin. 
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Fig. 5. — (a) Left panels: (top) IRAC 4.5 /im and 8.0 //m surface brightness profiles, averaged 
over elliptical contours; (bottom) color profile in units of [4.5]- [8.0] magnitudes. In the top 
panel, the 4.5 /im profile is shown as filled circles, and the 8.0 /im one as open circles. The 
exponential disk fit for R > 4" described in the text is shown as a solid line in the top panel, 
and the average disk color as a dotted line in the bottom one. The circumnuclear red "peak" 
at R ~ 4 " is an artifact due to the mismatch of the IRAC 4.5 and 8.0 /im point-response 
functions, (b) Right panels: The K-band surface brightness profile of Mrk996, averaged over 
elliptical contours (top panel). The three bottom panels show the J — H, H — K, and J — K 
color profiles. In these panels, the horizontal dotted line indicates the average color of the 
disk. The J — K color of the nuclear region is extremely red; most of this redness comes 
from the H — K color and less from the J — H color. 
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Fig. 6. — The 24 /zm (top left panel), 70 fim (top right), and 160 /zm (bottom) MIPS images 
of Mrk996 contoured on the 4.5/xm IRAC image. Contours run from the sky level + na 
to the maximum intensity, where n — 12, 5, and 3 for MIPS-24, MIPS- 70, and MIPS-160, 
respectively. The signature of the Airy ring, associated with a dominant point source, is 
evident in the 24 and 70 /im contours. 
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Fig. 7. — The IRS spectrum of Mrk 996 taken in the short- wavelength low- resolution mode 
and in both high-resolution modes: SL+SH spectra are shown in the upper panel, and the 
LH spectrum in the lower one. The PAH features at 5.7, 6.2, 7.7, 8.6, and 11.2 /im are clearly 
detected, as are several fine-structure emission lines: [Siv], [Nell], [Nelll], [S in], [Oiv]. 
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Fig. 8. — A blow-up of the IRS spectrum of Mrk996 around the [Oiv] 25.89 /xm line. [Fell] 
may be present, but is below the detection limit of the spectrum. 
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Fig. 9.— MMT spectrum of Mrk 996. 
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Fig. 10.— The (J - H,H - K) diagram. The colors of the nucleus of Mrk 996 and of 
two locations outside of it (filled squares from right to left connected by a dotted line) 
are compared to those of starburst galaxies (stars) (Hunt et al. 2002), blue compact dwarf 
galaxies (asterisks) (Hunt et al. 2002) and Seyfert galaxies (filled circles) (Alonso-Herrero 
et al. 2001). Mixing curves show how the colors change when various physical components 
contribute in increasing amounts to the observed disk emission (see Hunt & Giovanardi 1992). 
From top to bottom, the curves show the influence of extinction, of hot dust at 600 K, and 
at 1000 K, and of A stars; the end points of the mixing curves correspond to equal -fT-band 
contributions from stars and from the various components. The tick marks on the extinction 
line show magnitude increments in Ay units. 
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Fig. 11. — The SED of Mrk996 obtained by combining our UKIRT and Spitzer data with 
published data. The broadband data points are given by filled circles with error bars, upper 
limits by open circles, and IRS spectral points by x; the IRS spectrum has been rebinned 
into 1 [im. intervals. The global optical fluxes are taken from Thuan et al. (1996), assuming 
colors within a 9" aperture centered on the nucleus. Also shown are GRASIL (Silva et al. 
1998) models of Arp 220 (as a long-dashed curve) and M82 (dotted), scaled to the MIR flux 
between 20 and 30 /im. 



-39- 



t — i — rrrrra 




1 10 100 1000 

Rest-frame A (ycxm) 



Fig. 12.— Comparison of the Mrk996 SED with that of the Seyfert I galaxy NGC4151 
(Buchanan et al. 2006). 
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Fig. 13. — Fluxes F of the infrared lines (solid lines) listed in Table 1 and of some optical 
lines (dotted lines) computed with CLOUDY for a spherical high-density H II region with 
radius r. They are normalized to F max , the flux at log r = 21.25, where r is in cm. 
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Table I. IRAC and MIPS total fluxes of Mrk996 with Other Published Photometry 



Telescope / Instrument 


Wavelength 


Total flux 




(//m) 


(mJy) 


Spitzer/IRAC 


4.51 


5.8 ± 0.29 


» 


7.98 


12.7 ± 0.64 


Spitzer/MIPS 


23.7 


84 ± 4 




71.0 


381 ± 10 


;? 


156.0 


442 ± 12 


IRAS a 


12 


< 108 




25 


< 159 


)) 


60 


417 ± 50 


VI 


100 


< 713 


Uptical 


U.oo 


o i i n o 
Z.l ± O.Z 


» 


0.44 


3.7 ± 0.3 




0.55 


4.7 ± 0.4 




0.64 


6.9 ± 0.6 




0.79 


6.6 ± 0.6 


IRCAM/UKIRT C 


1.26 


9.6 ± 0.5 




1.65 


10.8 ± 0.5 




2.12 


8.9 ± 0.4 


2MASS a 


1.2 


10.2 ± 1 


» 


1.6 


12.6 ± 1 




2.2 


9.6 ± 1 



a From NED. 

b From Thuan et al. (1996), with colors defined 
within a 9" aperture centered on the nucleus. 

c In an aperture of 20" diameter. 



Table 2. Optical emission line intensities' 1 



Line 100xF(A)/F(H/3) 100x/(A)//(H/3) EW(H/3) 



3697 H17 


0.5 


± 


0.1 


0.7 


± 


0.1 


0.4 


3704 H16 


1.1 


± 


0.1 


1.6 


± 


0.1 


0.9 


3712 H15 


0.8 


± 


0.1 


1.2 


± 


0.3 


0.7 


3727 [0 ll] 


99.1 


± 


1.4 


144.4 


± 


1.8 


86.4 


3750 H12 


2.0 


± 


0.1 


2.9 


± 


0.3 


1.8 


3771 Hll 


2.5 


± 


0.1 


3.5 


± 


0.2 


2.2 


3797 H10 


3.3 


± 


0.1 


4.7 


± 


0.2 


2.9 


3820 He I 


1.7 


± 


0.1 


2.4 


± 


0.1 


1.5 


3835 H9 


5.3 


± 


0.1 


7.3 


± 


0.2 


4.7 


3868 [Ne ill] 


44.7 


± 


0.7 


61.7 


± 


0.8 


38.1 


3889 He I+H8 


13.8 


± 


0.2 


18.9 


± 


0.3 


11.9 


3968 [Ne lll]+H7 


24.7 


± 


0.4 


32.9 


± 


0.5 


21.4 


4009 He I 


0.4 


± 


0.1 


0.5 


± 


0.1 


0.3 


4026 He I 


2.4 


± 


0.1 


3.1 


± 


0.1 


2.0 


4038 N Il+N ll 


0.2 


± 


0.1 


0.2 


± 


0.1 


0.2 


4068 [S ll] 


1.4 


± 


0.1 


1.8 


± 


0.1 


1.2 


4076 [S n] 


0.3 


± 


0.1 


0.3 


± 


0.1 


0.2 


4102 m 


21.8 


± 


0.3 


27.8 


± 


0.4 


19.2 


4121 He I 


0.4 


± 


0.1 


0.4 


± 


0.1 


0.3 


4143 He I 


1.1 


± 


0.1 


1.4 


± 


0.1 


1.0 


4181 N ll 


0.3 


± 


0.1 


0.4 


± 


0.1 


0.3 


4207 ? 


0.7 


± 


0.1 


0.9 


± 


0.1 


0.6 


4232 I 


0.5 


± 


0.1 


0.6 


± 


0.1 


0.4 


4244 [Fe ll] 


0.7 


± 


0.1 


0.9 


± 


0.1 


0.7 


4267 C ll 


0.4 


± 


0.1 


0.6 


± 


0.1 


0.4 


4276 [Fe ll] 


0.5 


± 


0.1 


0.6 


± 


0.1 


0.5 


4287 [Fe ll] 


0.7 


± 


0.1 


0.9 


± 


0.1 


0.7 


4318 ll 


0.5 


± 


0.1 


0.6 


± 


0.1 


0.5 


4340 H 7 


42.7 


± 


0.6 


50.2 


± 


0.7 


40.6 


4363 [0 ill] 


14.3 


± 


0.2 


16.7 


± 


0.3 


13.6 


4388 He I 


0.6 


± 


0.1 


0.7 


± 


0.1 


0.6 



-43- 



Table 2 — Continued 



Line 100xF(A)/F(H/3) 100x/(A)/7(H/3) EW(H/3) 



4415 


[Fe nj 


1.0 


± 


0.1 


1.2 


± 


0.1 


1.0 


4452 


[Fe II J 


0.5 


± 


0.1 


0.6 


± 


0.1 


0.5 


4471 


He I 


6.1 


± 


0.1 


6.8 


± 


0.1 


5.9 


4520 


Si II 


0.4 


± 


0.1 


0.5 


± 


0.1 


0.4 


4658 


[Fe ill] 


2.8 


± 


0.1 


3.0 


± 


0.1 


2.6 


4702 


[Fe ill] 


0.6 


± 


0.1 


0.7 


± 


0.1 


0.6 


4713 


[Ar lv]+He I 


1.2 


± 


0.1 


1.2 


± 


0.1 


1.1 


4734 


[Fe ii] 


0.6 


± 


0.1 


0.6 


± 


0.1 


0.6 


4755 


[Fe ill] 


0.5 


± 


0.1 


0.5 


± 


0.1 


0.5 


4770 


[Fe ill] 


0.2 


± 


0.1 


0.2 


± 


0.1 


0.2 


4778 


[Fe ill] 


0.3 


± 


0.1 


0.3 


± 


0.1 


0.3 


4815 


[Fe ii] 


0.4 


± 


0.1 


0.4 


± 


0.1 


0.4 


4861 


H/3 


100.0 


± 


1.5 


100.0 


± 


1.5 


106.9 


4881 


[Fe ii] 


0.4 


± 


0.1 


0.4 


± 


0.1 


0.4 


4903 


[Fe ii]+Si II 


2.0 


± 


0.1 


2.0 


± 


0.1 


2.2 


4922 


He I 


2.2 


± 


0.1 


2.2 


± 


0.1 


2.4 


4959 


[0 m] 


100.9 


± 


1.5 


98.2 


± 


1.4 


110.3 


4987 


[Fe ill] 


0.5 


± 


0.1 


0.5 


± 


0.1 


0.6 


5007 


[0 m] 


300.0 


± 


4.3 


287.9 


± 


4.2 


328.4 



a C(H/3) = 0.53, F(H/3) = 7.5xl0- 17 W m" 2 , J(H/3) = 25.4xl0- 17 W 
in" 2 , EW abs = A. 
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Table 3. PAH Parameters Obtained by Fitting Lorentzian Profiles 



Wavelength a 


Integrated 15 


Equivalent 13 


FWHM b 




Flux 


Width 




(H 


(10- 17 Wm- 2 ) 


(/iHl) 


(/iHl) 


5.759 


3.3 (0.1) 


0.245 (0.029) 


0.044 (0.002) 


6.225 


14.8 (0.4) 


1.030 (0.046) 


0.172 (0.005) 


7.771 


54.8 (2.6) 


3.640 (0.427) 


0.837 (0.025) 


8.664 


6.1 (0.4) 


0.415 (0.064) 


0.200 (0.008) 


11.278 


21.9 (1.0) 


1.360 (1.011) 


0.450 (0.260) 



a Rest wavelength, corrected for z = 0.00541. 

b Standard deviations of the repeated measurements are 
given in parentheses. The true uncertainty including calibra- 
tion is probably ~20%. 



-45 - 



Table 4. Fine Structure Line Parameters Obtained by Fitting Gaussian Profiles 



Line a 


E- b 

- LJ ion 


Nominal 


Fitted c 


Integrated* 3 


Equivalent 


FWHM d 






Wavelength 


Wavelength 


Flux 


Width 






(eV) 




(/iin) 


(10- 17 Wm- 2 ) 


(/iin) 


(H 


[Siv] 


34.8 


10.511 


10.512 


3.4 (0.15) 


0.27 (0.14) 


0.026 (0.0004) 


[Neil] 


21.6 


12.814 


12.814 


3.8 (0.3) 


0.14 (0.04) 


0.082 (0.062) 


[Nem] 


41.0 


15.555 


15.556 


10.77 (0.5) 


0.22 (0.01) 


0.035 (0.005) 


[Sin] 


23.3 


18.713 


18.714 


7.4 (0.5) 


0.13 (0.01) 


0.039 (0.0006) 


[Oiv] 


54.9 


25.890 


25.901 


1.2 (0.3) 


0.02 (0.01) 


0.026 (0.005) 


[S in] 


23.3 


33.481 


33.480 


10.85 (4.8) 


0.24 (0.10) 


0.055 (0.035) 



a Because [Ar ill] and [Ar il] are seen only in the low-resolution spectra, we do not give their 
parameters here. 

b Ionization potential. 

°Rest wavelength, corrected for z = 0.00541. 

d Standard deviations of repeated measurements are given in parentheses. The true uncer- 
tainty including calibration is probably ~20%. 



Table 5. Parameters of narrow and broad emission lines 3. 



Narrow component Broad component 



Line 


100xF(A)/F(H/3) 


FWHM (A) 


100xF(A)/F(H/3) 


FWHM (A) 


3868 [Ne ill] 


38.1 


3.5 


52.5 


8.5 


3889 He I+H8 


16.6 


2.6 


11.8 


6.1 


4102 m 


23.8 


2.6 


20.7 


6.7 


4340 H 7 


42.2 


2.6 


42.5 


6.7 


4363 [0 ill] 






31.2 


7.6 


4471 He I 


4.3 


2.8 


8.6 


8.4 


4686 He II 






11.4 


31.7 


4861 H/3 


100.0 


2.6 


100.0 


7.6 


4959 [0 ill] 


126.0 


2.7 


69.2 


6.6 


5007 [O in] 


399.6 


2.9 


175.9 


7.5 



a F(R(3) nar = 4.0xl0~ 17 W m" 2 , F(R(3) br = 3.5xlQ- 17 W m" 2 . 
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Table 6. Comparison of observed and modeled optical emission line fluxes 3. 



Line 


Observ. 


CLOUDY Models b 
Model I Model Ia c 


Model I/Obs. 


Model Ia/Obs. 


3727 [0 II] 


36.7 


73.3 


33.8 


2.00 


0.92 


3868 [Ne ill] 


15.7 


15.0 


14.5 


0.95 


0.93 


4363 [0 ill] 


4.3 


5.1 


5.1 


1.18 


1.18 


4861 H/3 


25.4 


41.8 


28.4 


1.65 


1.11 


5007 [0 ill] 


73.1 


64.7 


64.2 


0.89 


0.88 



a Fluxes are in units of 10~ 17 W m~ 2 . 

b Fluxes are calculated assuming the distance D — 21.7 Mpc. 

c Model la is the same as Model I except that the radius of the modeled H II 
region is 410 pc (log r = 21.1) instead of 580 pc. 
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Table 7. Comparison of observed and modeled MIR emission line fluxes 3. 



CLOUDY Models b Model/Observ. 



Line 


Obs. 


Mod.I 


Mod.II 


AGN 


WNE-w 


Shock c 


A d 


B e 


C f 


Ds 


10.5 [S iv] 


3.42 


0.19 


2.49 


0.65 


0.80 


0.62 


0.78 


0.97 


1.02 


0.96 


12.8 [Ne ll] 


3.81 


4.21 


0.54 


0.39 


0.03 


0.05 


1.25 


1.35 


1.25 


1.26 


15.6 [Ne ill] 


10.77 


3.32 


6.88 


1.08 


1.02 


0.88 


0.95 


1.05 


1.04 


1.03 


18.7 [S in] 


7.39 


2.90 


3.66 


0.66 


0.38 


0.39 


0.89 


0.98 


0.94 


0.94 


24.3 [Ne v] 










0.34 


0.12 


0.18 










25.9 [0 iv] 


1.19 








1.12 


1.05 


0.93 





0.94 


0.88 


0.78 


33.5 [S in] 


10.85 


3.33 


5.21 


0.88 


0.55 


0.56 


0.79 


0.87 


0.84 


0.84 



a Fluxes are in units of 10~ 17 W m~ 2 . 

b Fluxes are calculated assuming the distance D = 21.7 Mpc. 

c Shock velocity is 250 km s -1 . 

d A=I+II. 

C B=I+II+AGN. 

f C=I+II+WNE-w. 

g D=I+II+Shock. 



